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^s^j ■ Abstract. A small fraction of all quasars are strongly lensed and multiply imaged, with usually a galaxy acting 

' as the main lens. Some, maybe all of these quasars are also affected by microlensing, the effects of stellar mass 

objects in the lensing galaxy. Stellar microlensing not only has photometric effects (the apparent magnitudes of 
the quasar images vary independently due to the relative motion between source, lens and observer) , it also affects 
the observed position of the images. This astrometric effect was first explored by Lewis and Ibata (1998): the 
position of the quasar - i.e. the center-of-light of the many microimages - can shift by tens of microarcseconds 
due to the relatively sudden (dis-)appearance of a pair of microimages when a caustic is being crossed. 
We explore this effect quantitatively for different values of the lensing parameters n and 7 (surface mass density 
and external shear) covering most of the known multiple quasar systems. We show examples of microlens-induced 
' quasar motion and the corresponding light curves for different quasar sizes. We evaluate statistically the occurrence 

of large shifts in angular position and their correlation with apparent brightness fluctuations. We also show 
statistical relations between positional offsets and time from random starting points. As the amplitude of the 
astrometric offset depends on the source size, astrometric microlensing signatures of quasars - combined with the 
photometric variations - will provide very good constraints on the size of quasars as a function of wavelength. 
We predict that such signatures will be detectable for realistic microlensing scenarios with near future technology 
in the infrared/optical (Keck-Interferometry, VLTI, SIM, GAIA). Such detections will show that not even high 
CZ3 . redshift quasars define a "fixed" coordinate system. 

Key words, cosmology - gravitational lensing - quasars - astrometry 

■ 1. Introduction Lewis and Ibata (1998): when the quasar crosses a (micro- 
caustic, a new pair of microimages is created or destroyed. 

Gravitational lensing acts on quasars in a num- The relatively sudden (dis-) appearance of a highly magni- 

ber of ways. The best known effect is strong lens- fled image pair does not just produce a strong fluctuation 

ing which produces multiply imaged quasars. Only a in apparent brightness, it can also shift the center-of-light 

small fraction of all quasars (about 1 in 500) arc of the quasar (the weighted sum of all the microimages) 

multiply imaged with galaxies acting as the main by tens of microarcseconds. Such a positional change is in 

lense s (for an updated list, see the CASTLES web principle detectable with current VLBI techniques in the 

page http://cfa-www.harvard.edu/glensdata/). Most radio regime, and will be measurable in the infrared and 

of these quasars - and possibly some "single" quasars optical wavebands as well with (near) future technology 

as well - are also affected by microlensing: the coher- (Keck-Interferometry, VLTI, SIM, GAIA). 
ent effects of stellar mass objects in the lensing galaxy. 

Microlensing is well known for affecting the apparent A number of studies have explored astrometric mi- 
brightness of the quasar, which changes due to the rel- crolensing in a different regime, namely stars in the Milky 
ative motion between lens, source and observer (see, Way (or dark matter objects in its halo) acting on back- 
e.g., Paczyhski 1986, Wambsganss, Paczyhski & Schneider ground stars either in the Magellanic Clouds or in the 
1990, Lewis et al. 1998, Wyithe et al. 2000a, Wambsganss bulge of the Milky Way (e.g., Miyamoto & Yoshii 1995; 
2001, Wozniak et al. 2000a,b). However, stellar microlens- Miralda-Escude 1996; Mao & Witt 1998; Boden, Shao & 
ing has yet another effect, which was first explored by van Buren 1998; Goldberg & Wozniak 1998; Paczyhski 



2 



Marie Treyer and Joachim Wambsganss: Astrometric Microlensing of Quasars 



1998; Han, Chun & Chang 1999; Han & Kim 1999; Han 
& Jeong 1999; Safizadeh, Dalai, Griest 1999; Dominik & 
Sahu 2000; Gould & Han 2000; Han & Kim 2000; Salim & 
Gould 2000; Delplancke et al. 2001; Belokurov & Evans 
2002; Dalai & Lane 2003). The effects of Milky Way 
stars on background quasars were studied by Hosokawa, 
Ohnishi & Fukushima (1997), Sazhin et al. (1998) and 
Honma & Kurayama (2002). Williams & Saha (1995) had 
discussed large image shifts produced by substructure in 
the lensing galaxy. However, not much work has yet been 
done on cosmological astrometric microlensing - stars in 
lensing galaxies acting on even more distant quasars - be- 
yond Lewis & Ibata (1998), with the recent exception of 
Salata & Zhdanov (2003). The present paper aims to ex- 
plore this further. 

We study quantitatively eight different cases with 
various values of the surface mass density, with and 
without external shear. We present example microlens- 
ing lightcurves with the corresponding center-of-light 
shifts and we investigate the correlation between high- 
magnification photometric events and large-offset astro- 
metric events. 

As first pointed out by Lewis & Ibata (1998), large off- 
set events are typically correlated with high-magnification 
events, whereas the inverse is much less true. The reason 
for this is that the location of the newly appearing bright 
image pair during a caustic crossing is unrelated to the 
previous center-of-light of the microimages. In rare cases, 
the new image pair may appear at or close to the center- 
of-light of the pre-existing microimages. Such a situation 
would correspond to a large change in brightness with very 
little change in position. In most cases, however, the new 
bright image pair will appear at a location which is unre- 
lated to the previous center-of-light, and hence produce a 
sudden jump. Since the positional offset is preferentially 
perpendicular to the direction of the external shear, the 
shear direction can be inferred this way. 

After introducing the microlensing length and time 
scales (Section 2), we describe our simulations (Section 
3). In Section 4, we illustrate the effect of astrometric and 
photometric fluctuations and present statistical correla- 
tions between positional offsets, magnitude fluctuations 
and time intervals between the measurements for four mi- 
crolensing situations with different surface mass densities 
(k = 0.2,0.4,0.6,0.8), with and without external shear 
(7 = k or 0). Finally we discuss the possibilities of real 
detections of this phenomenon in the near future. 

2. Microlensing basics: length, time and angular 
scales 

2.1. Standard mass, length and time scales 

The lensing effects of cosmologically distant compact ob- 
jects in the mass range 10~ 3 < M/M < 10 3 on back- 
ground sources is usually called "cosmological microlens- 
ing". The source is typically a distant quasar, but in 
principle other objects can be microlensed as well, e.g. 



high-rcdshift supernovae (e.g., Rauch 1991) or gamma- 
ray sources/bursters (Torres et al. 2003; Koopmans & 
Wambsganss 2001). The only condition is that the source 
size be comparable to or smaller than the Einstein radius 
of the foreground lens. 

The microlenses can be ordinary stars, brown dwarfs, 
planets, black holes, molecular clouds, or other compact 
mass concentrations. In most practical cases, the mi- 
crolenses are part of a galaxy which acts as the main 
(macro-) lens. However, microlenses could also be located 
in, say, clusters of galaxies (Tadros, Warren & Hewett 
1998, Totani 2003) or they could even be imagined "float- 
ing" freely and filling intergalactic space (Hawkins 1996, 
Hawkins & Taylor 1997). 

The relevant length scale for microlensing is the 
Einstein radius in the source plane: 
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where "typical" lens and source redshifts of zl = 0.5 and 
zs — 2.0 were assumed for the expression on the right 
hand side (G and c are the gravitational constant and the 
velocity of light, respectively; M is the mass of the lens, 
Dl, Ds, and Dls are the angular diameter distances be- 
tween observer - lens, observer - source, and lens - source, 
respectively, and a concordance cosmological model is as- 
sumed With fi tot = 1, ^matter = 0.3, £1 A = 0.7). 

This length scale translates into an angular scale of: 
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It is obvious that image splittings on such small angu- 
lar scales can not be observed directly. What makes mi- 
crolensing observable in the first place is the fact that 
observer, lens(es) and source move relative to each other. 
Due to this relative motion, the microimage configuration 
changes with time, and so does the total magnification, i.e. 
the combined fluxes of all the microimages which make up 
the macro-image. This change in magnification over time 
(lightcurve) can be measured: microlensing is a dynamical 
phenomenon. 

There are two time scales involved: the standard lens- 
ing time scale ts is the time it takes the lens (or the source) 
to cross a length equivalent to the Einstein radius, i.e. 
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where the same typical assumptions are made as above, 
and the relative transverse velocity 1 t>6oo is parametrized 
in units of v± = 600km/sec and for a Hubble constant 
Hq = 75/175 km/sec/Mpc. This time scale results in rela- 
tively large values. However, microlensing fluctuations are 



1 For simplicity we assume here that all the transverse mo- 
tion is done by the lensing galaxy 
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Table 1. Summary of scales for the "typical" lensing case [z L ~ 0.5, zq k, 2) and for the special case of Q2237+0305 
(z L w 0.039, zq ~ 1.69): 

typical lensed quasar case lensed quasar Q2237+0305 

Physical Einstein radius in quasar plane R E 3.2 x 10 16 y/M/M® hf 5 a5 cm 1.6 x 10 16 yj 'M/M Q h^ 0,5 cm 

Angular Einstein radius 6 E 2.2yJ M/M Q hf 5 ' 5 /iarcsec 7.3y/ M/M® hf 5 ' 5 parcsec 

Einstein time t E 25y/M/M® v^ Q hf 5 as years 8.7y/M/M® v^ hf 5 as years 

Crossing time t crosa 210 Ris v$q 1% ' 5 days 18 R15 v^q h^g ' 5 days 



expected (and observed, see Wozniak et al. 2000a, b) on 
much shorter time scales. The magnification distribution 
is highly non-linear with sharp caustic lines separating re- 
gions of low and high magnification. Often, the density 
of caustics is quite high, so that a source may encounter 
half a dozen or more caustic lines within the length of 
one Einstein radius. When a source crosses a caustic, we 
can expect a large change in magnification (and corre- 
spondingly in position) within the time it takes the source 
to cross its own radius. So the relevant time scale is the 
"crossing time" : 



, /l I \ Rsource 

Across 1,1 + Z L) r 
v±(D s /D L ) 



0.58 R 15 u 6o6 h 75 5 years 
210 R 15 v m l h7 5 5 days. 



(4) 



Here the quasar size i?is is parametrized in units of 
10 15 cm. 



2.2. The special case of Q2237+0305 

The quadruple quasar Q2237+0305 (Huchra et al. 1985; 
Irwin et al. 1989; Wambsganss et al. 1990; Wyithe et al. 
2000a, b; Wozniak et al 2000a, b) is a very special and 
favorable case and of particular interest to microlensing 
studies. It was the first system in which microlensing was 
discovered (Irwin et al. 1989). Subsequently it received a 
lot of attention, both observational (Corrigan et al. 1991, 
Ostcnsen et al. 1996; Wozniak et al 2000a, b) and theo- 
retical (Wambsganss et al. 1990; Wyithe et al. 2000a, b; 
Yonehara 2001). Due to the fact that the lensing galaxy 
is so close (zq — 0.039, Huchra et al. 1985), the physical 
and angular Einstein radii are considerably different from 
the standard case treated above: 
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6^2237+0305 « 7.3 X 10 V ^ h 75°' 5 arCSeC > ( 6 ) 

The resulting time scales (Einstein time and crossing 
time) are much shorter than in almost all other multiple 
quasars: 
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For that reason, this quadruple system is ideally suited for 
microlensing studies. The length, time and angular scales 
for the "typical case" as well as for the special case of 
Q2237+0305 are summarized in Tabled 

3. The Simulations 

In order to explore astrometric microlensing for a vari- 
ety of realistic scenarios, we consider eight different cases, 
with the following values of the dimensionless surface mass 
density: k = 0.2, 0.4, 0.6 and 0.8. Each one of these we 
treated both without external shear (7 = 0) and with ex- 
ternal shear equal to the surface mass density (7 = k), 
corresponding to an isothermal sphere model for the lens- 
ing galaxy. The shear always acts along the X-axis in the 
simulations. This means that the caustics are expanded 
in the X-direction and compressed in the Y-direction, as 
can be seen for example in the top left panel of Fig. 
Horizontal bands appear on large scales (many Einstein 
radii). The effect is not always obvious on small scales 
however (sub-Einstein radii, bottom left panel of Fig.[3J). 
Inversely, if we follow the light rays forward from source to 
observer (as opposed to backwards as in the simulations, 
which is explained below), the micro-image configuration 
of a source will appear compressed in the X-direction and 
expanded in the Y-direction (cf. the right hand panels of 
Fig. El. 

For each of these eight cases, we produced a magnifi- 
cation pattern in the source plane with a side length of 
La = 20 Re, sampled on a 1000 2 pixel grid, i.e. 1 Re is 
covered by 50 pixels, or 1 pix = 0.02 Re- This is our 
starting configuration A. The minimum source size we 
can consider is given by the pixel scale. In order to ex- 
plore smaller source sizes as well (at the cost of smaller 
magnification patterns which may not be representative), 
we zoom in into the central part of the first configura- 
tion. We do higher resolution simulations (all on 1000 2 
pixel grids) in three steps of factors of two (configurations 
B, C and D, respectively), resulting in magnification pat- 
terns with side lengths of Lb = IORe, Lc — 5Re, and 
Le = 2.5Re, respectively. For the highest resolution sim- 
ulation (configuration D), this means Re — 400 pixels, or 
1 pix = 0.0025i? B . 

We used a modified version of the ray shooting code 
described in Wambsganss (1990, 1999). In the original 
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Fig. 1. Illustration of source position in magnification 
pattern/source plane (left) and corresponding microim- 
age configuration in the image plane, for cases k = 0.2, 
k = 0.4, k = 0.6, and k = 0.8 (from top to bottom) and 
no external shear (7 = 0.0). The center-of- light position 
is marked with a plus sign. The horizonal bar indicates a 
length of two Einstein radii (2 Re)- Note the change in 
scale in the right hand panels. 

version, rays are followed backwards from the observer 
through the lens plane (where all the deflections are de- 
termined) to the source plane (where the light rays are 
"collected" in pixels). The density of rays then indi- 
cates the magnification as a function of position in the 
source plane, often displayed as color coded magnifica- 
tion patterns. Lightcurves can be obtained by convolving 




Fig. 2. Same as Figure 1, here for the cases with external 
shear in the X-axis direction: n = 7 = 0.2, = 0.4, = 0.6, 
and = 0.8 (from top to bottom); the center-of-light posi- 
tion is marked with a plus sign, the horizonal bar indicates 
a length of two Einstein radii. 



a given source profile with this two-dimensional magnifica- 
tion map. However, all information about where the rays 
originated from in the lens/image plane is lost in this orig- 
inal algorithm. For the exploration of the centroid shift of 
the collection of microimages, it is exactly this information 
that is required. Hence we modified the code to record the 
positions of all the individual microimages brighter than 
a given magnification threshold. To do this, we defined a 
second, larger regular grid of 4000 2 rays, which covers in 
the image plane a region of four times the angular side 
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Table 2. Source sizes (Gaussian widths) used in our mi- 
crolensing simulations (configurations A and B) in units of 
pixels, Einstein radii and in physical units for the typical 
lensed quasar and for the special case of Q2237+0305 
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(typical case) 


(Q2237+0305) 
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16 


0.32 


1.2 x 10 16 


4 x 10 1B 
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0.16 
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2 x 10 16 
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1 x 10 16 
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5 x 10 15 
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0.04 


1.5 x 10 15 


5 x 10 15 
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0.02 


7.5 x 10 14 


2.5 x 10 15 



length of the magnification pattern (i.e. (80 9e) 2 for con- 
figuration A). For a set of 1.6 x 10 7 test rays, we keep 
track of the positions in the image plane as well as in the 
source plane. In this way, we can find all the positions in 
the image plane that are mapped onto a certain area in 
the source plane, and hence identify all the microimages 
corresponding to a particular source position and size. 

In order to get information on both the magnification 
and the microimage locations of a finite source at a given 
position, we convolved the magnification pattern with a 
luminosity profile. We used circular sources with Gaussian 
widths er = 2, 4, 8 and 16 pixels, corresponding to phys- 
ical sizes between 0.04 Re and 0.32 Re in the starting 
configuration A, or 1.5 x 10 15 cm to 1.2 x 10 16 cm in the 
typical case described in Section 2.1. We did the same 
with the higher resolution magnification patterns, config- 
urations B, C and D. However, it turned out that the high 
resolution cases, although well suited for studying individ- 
ual caustic crossings for small sources, are not quite large 
enough for statistical investigations. For this reason, we 
restricted ourselves to cases A and B for the statistical 
evaluations below. The numerical values corresponding to 
the different source sizes we used are tabulated in Tabled 

Using these simulations, we determined the positions 
of the individual microimages, the center-of- light and the 
total magnification of the macro-image corresponding to 
a particular source position and profile. In Figs.HandGl a 
microimage situation is shown for each of the eight cases 
considered: k = 0.2,0.4,0.6 and 0.8 with 7 = (Fig.UJ), 
and with 7 = k acting along the X-axis direction (Fig. 
[21). In the left columns, a small part of the magnification 
pattern is shown, with the source position and profile su- 
perimposed. The panels on the right hand side show the 
particular microimage configurations, with the light cen- 
troid indicated by a plus sign (cf. Paczynski 1986). 

In Figure |3 the method is illustrated for n = 0.8 
and 7 = 0.8. The top left panel shows the magnifica- 
tion pattern in the source (quasar) plane with side length 
L = 10 6e (B configuration). The grey scale indicates 
regions of different magnifications in the source plane: 
the lighter the grey, the higher magnification. The large 
square indicates the region in the source plane for which 




Fig. 3. Illustration of the ray shooting method, here for 
case k = 0.8, 7 = 0.8. Top left: Magnification pattern in 
the source (quasar) plane with side length L — 10 Be (B 
configuration). The large square indicates the region in the 
source plane for which all source positions were evaluated. 
Top right: Corresponding microimage configuration of the 
large square "source" on the left, side length is 40 9e- 
These are the parts of the sky (image plane) where light 
bundles originated from the square. Note the strong effect 
of shear in the Y-axis direction. Bottom left: The square 
here corresponds to the small square in the top panel, 
zoomed in 8 times. Bottom right: Corresponding micro- 
image configuration for the square on the left. 



all source positions were evaluated. The top right panel 
shows the corresponding microimage configuration in the 
image plane: the white regions are the parts of the image 
plane where light bundles appear which originated from 
within the square on the left (side length is 40 9e)- In 
other words, it shows what a large square shaped source 
would look like to the observer. The bottom panels show 
the same for the smaller square region indicated in the top 
left panel, zoomed eight times. The many isolated light 
patches indicate that microimages are spread over a very 
large area in the image plane. 

Figure 0] shows a quasar microlensing scenario with 
microlensing parameters k = 0.6 and 7 = 0.6, and side 
length 10 8e (B configuration). The straight vertical white 
line marks the track of the quasar motion relative to the 
magnification pattern; the length of the path is 2.0 Re- 
For this particular track (followed from the lowest part up- 
wards), Figure [S] shows from top to bottom: the X- and Y- 
coordinates of the quasar relative to the starting position 
(A9x and A9y) as a function of time; the absolute value of 
the positional shift ((A^l) relative to the starting position 
as a function of time; and the corresponding light curve 
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Fig. 4. Illustration of a quasar microlensing scenario (here 
for lensing parameters n = 0.6 and 7 = 0.6, sidelength 
about 10 Re)- The gray scale indicates regions of different 
magnifications in the source plane. The straight vertical 
line marks the track of the quasar relative to the magnifi- 
cation pattern. The total length of the track is 2.0 Re- 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

t/t E 



Fig. 5. The four lines indicate from top to bottom: a) 
A8x , the X-coordinate of the quasar relative to the start- 
ing position as a function of time; b) same for the Y- 
coordinate; c) the absolute value of the positional shift 
|A0| relative to the starting position as a function of 
time; and d) the corresponding light curve of the quasar 
for the example track shown in Figures 0] The solid 
lines correspond to a sourc e size a = 0.02Re (typically 
about 6.4 x 1 14 y/M /M Q ' 5 cm) the dashed lines to 
a = O.WRe (typically about 5.1 x 10 1 5 ^M/M Q h^- 5 cm). 



(Am) of the quasar. The solid and dotted lines correspond 
to two different values of the source size: a — O.QARe (4 
pixels in B configuration) and a = 0.16i?E (16 pixels in B 
configuration) , respectively. The track in Fig. 0] starts in 
a region of low magnification which is taken as the zero 
point of the magnitude scale on the lowest panel in[3 

The two panels in Figure HJ1 represent the centroid shift 
for the two different source sizes. The whole track has a 
length of 2.0 tg (^50 and 15 years for the 'typical' lensing 
case and the Q2237+0305 case, respectively). The labels 
to, t\ and ti correspond to the starting, middle and fi- 
nal positions, respectively. Snapshots of the ensemble of 
microimages at times to, t\ and ti are shown in Figure 
assuming a = Q.16Re- All microimages are plotted up 
to a radius of 3 a, without accounting for the declining 
source profile. On average, the brightness of the macroim- 
age declines as the fourth power of distance to the center 
of light (Katz, Balbus, Paczyfiski 1986). On our plots, due 
to the finite resolution, some of the distant and faint mi- 
croimages appear bigger than they are. This figure simply 
aims at illustrating the very large spatial spread of the 
microimages, covering roughly 30 Be X 40 Be- The lower 
panels show enlargements of the central parts with the 
path of the center-of-light superimposed. 



4. Results and Discussion 

4.1. Microlens-induced positional offsets and 
corresponding magnitude changes 

In order to evaluate the correlations between the magni- 
tude changes and the positional changes of a microlensed 
quasar, we simulated pairs of measurements separated 
by time intervals of 0.02 t E (Figures B El and 0.16 t E 
(Figures HID , respectively. These values correspond to 
about half a year and four years in the "typical" case, and 
to two months and 1.5 years in the case of Q2237+0305 
(cf. Table QJ. We placed the source at a random position 
in the magnification pattern and chose a second position 
at a distance of 0.02J?e or 0.16i?E either parallel or per- 
pendicular to the action of the external shear (X-axis) . We 
determined differences in the magnifications |Am| and the 
center-of-light positions | AB\ between these two source po- 
sitions. Each pair of measurements (|Ato|, |A0|) is repre- 
sented as a point in the various panels of Figures|Slto^2 In 
each panel, the offset in position |A0| is shown against the 
corresponding offset in magnitude |Am| for about 20,000 
such measurement pairs. In the left columns of all four 
figures, a small Gaussian source size of a — 0.02Re is as- 
sumed, in the right column the source size is a = 0.16Re- 
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Fig. 6. Effect of the source size on the centroid shift for 
the example track shown in Fig. 01 and (k = 0.6 and 
7 = 0.6). The two panels show the astrometric tracks that 
would occur for the linear quasar motion shown in Fig. 0] 
for two different source sizes: a = 0.04 and 0.16 Re (4 
and 16 pixels in configuration B, respectively). The total 
length of the track is 2.0 or about 15 years for the 
quasar Q2237+0305. The labels to, t\ and t 2 correspond 
to the starting, central and final positions, respectively. 



In each panel, the thick (lower) line indicates the median 
of |A0| as a function of |Am|, and the thin (upper) line 
shows the 95th-percentile: 5% of all simulations for a given 
| Am | would result in offsets | A8\ which are above this thin 
line. 

In Figure [SI the four lensing situations without exter- 
nal shear are considered: k = 0.2,0.4,0.6,0.8 (from top 
to bottom) with a small time step At — 0.02^. In all 
panels, the majority of the points cluster near the origin. 
This is easily understandable: the position of the source 
relative to the caustics has not changed by much during 
the short time interval, so the changes both in magnifi- 
cation and in the centroid position tend to be small. For 
small sources (left panels), however, |Am| and |A0| can 
reach relatively large values with the median lines indi- 
cating an almost linear statistical relation. The slope of 
these median lines slightly increases with increasing sur- 
face mass density k. The time step corresponds roughly to 
the crossing time: At = 0.02^ w t cross , and there are in- 
deed cases with easily measurable magnitude fluctuations 
(|Am| > 0.2 mag), which result in center-of-light offsets 
of about 0.5 6e- In the right panel, there are no significant 
changes in either magnification or center-of-light position, 
because the source has moved only a fraction of its own 
diameter (At < t CIOSS ). 



Fig. 7. Top: Three snapshots of the ensemble of microim- 
ages assuming er = 0.16i?£. All bright microimages are 
plotted up to a radius of 3 a (due to the finite resolu- 
tion of the plotting routine, some of the faint microim- 
ages are represented larger than they actually are. Left 
to right: starting, middle et final positions of the track in 
Fig- IH corresponding to t=0.0, 1.0 and 2.0 tE, respectively 
(note the angular scale: the microimage distribution covers 
about 30 9e x 40 6e\)- Bottom: zoomed central regions; 
the respective positions of the center-of-light for the three 
epochs are marked as crosses on the track of the center- 
of-light (with the very central source parts corresponding 
to 0.5 a indicated in red). 

In Figure the same is shown for a larger time step 
At = O.lGtE- Many more points are now spread towards 
larger offsets and larger magnification changes. For the 
small source, median values of |A0| > We are reached in 
the highest surface mass density cases. In fact, the 95th- 
percentile line for the k = 0.8 case indicates that for mag- 
nitude changes |Am| > 0.4mag, 5% of the all cases re- 
sult in center-of-light offsets larger than 49 e- For the large 
source (right column) the expected offsets are still quite 
moderate, with median values of about 0.5 Be- This is not 
too surprising, because the time interval corresponds to 
just about the crossing time for the source. 

Figures EH an d E] contain the same diagrams for the 
cases with external shear, n = 7 = 0.2,0.4,0.6,0.8 (from 
top to bottom). Larger offsets are reached here than in 
the corresponding scenarios without shear. Especially the 
two middle rows with k — 7 = 0.4 and 0.6 - which best 
represent the typical values of convergence and shear of 
a multiply imaged quasar - produce values of positional 
offset a factor of two higher than the corresponding cases 
without shear: For short time steps and small sources, 
medians of |A0| « We and 95th-perccntiles of |A0| sa 
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Fig. 8. Offset in position |A0| versus offset in apparent 
brightness |Am|, for two simulated measurements sepa- 
rated by At — 0.02tE- This time interval corresponds 
to about half a year for the "typical" lensing case and 
to two months for the case of Q2237+0305 (cf. Table 0. 
From top to bottom, the panels represent the cases n = 
0.2, 0.4, 0.6 and 0.8, respectively, without external shear 
(7 = 0.0). The source sizes are a — Q.02Re (left column), 
and O.WRe (right column). The median of the |A0| dis- 
tribution as a function of |Am| is shown as a thick line 
(points marked as crosses), the 95th-percentile is indicated 
as a thin line (and squares) . 

28 e (small source, left column) are reached. For smaller 
values of the surface mass density, the caustic density and 
hence the number of microimages is not large enough to 
produce big positional offsets, whereas for higher values 
of k, the density of caustics is so high that an additional 
microimage pair only produces a small fluctuation. 

Results are even more dramatic when assuming a 
larger time step (At = 0.16ts), as displayed in Figure 
1111 median values of |A0| between 0.58 e and We for the 
small source size (left column), and significant median val- 
ues of |A0| « 0.58e to 19e even for the large source (right 
column) . 

4.2. Microlens-induced positional offset as a function 
of time (with magnitude thresholds) 

We evaluated the shifts in center-of-light positions for 
time intervals intermediate between 0.02 < At/ts < 0.16 
assuming a threshold for the magnitude fluctuations. In 
Figures IT51 and IT51 the median and 95th-percentile offsets 
are shown as a function of increasing time step At, for 
I Am I > 0.2 mag (left panels) and |Am| > 0.5 mag (right 
panels). This was done with the following idea in mind: 




Fig. 9. Same as Fig. |SJ but here for two measurements 
separated by a larger time interval At = O.lQtE (corre- 
sponding to about four years in the "typical" case, and 
1.5 years for the case of Q2237+0305, cf. Tabled]): source 
sizes are <j = 0.02Re (left column), and 0.16i?E (right 
column); median of |A0| as function of |Am| is thick line 
(points are crosses), 95th-percentile is thin line (squares). 



ic^O.2, '7-0.2. ij^O.16 R t . ai-0.02 1, 




*='9,i.V'S©.<, o O.'.fi S„, 41=0.02 l. n 
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' 1 I ' ' ' I 1 1 1 L 1 1 1 I ' ' 

«-0.B, 7=0.8. o=O.I6 R r . AL 0.02 l. c 



02 04 0.6 D.B 1 02 0.4 0.6 0.8 1 

Amag Am ag 

Fig. 10. Same as Fig. [5J but here for the cases with ex- 
ternal shear 7 = k = 0.2,0.4.0.6,0.8 from top to bottom: 
source sizes are a — Q.02Re (left column), and 0.16-Re 
(right column); median of |A0| as function of \Am\ is 
thick line (points are crosses), 95th-percentile is thin line 
(squares). 
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Fig. 11. Same as Fig. El but here for the two measure- 
ments separated by a larger time interval At = 0.16<e: 
source sizes are a — 0.02Re (left column), and 0.16i?£ 
(right column); median of |A0| as function of |Am| is 
thick line (points are crosses), 95th-percentile is thin line 
(squares) . 



since very accurate measurements of quasar image posi- 
tions are "expensive" , it is unlikely that all lensed quasar 
candidates can be astrometrically monitored. However, 
photometric monitoring is comparably cheaper. So ideally, 
one could determine the positions of the most promising 
multiple quasars with high accuracy once, then monitor 
them photometrically, and whenever a large microlens- 
induced magnitude change has been detected, a second 
astrometric measurement should be performed. 

In each panel of Figure ^1 the two sets of curves show 
the median (thick lines) and the 95th-percentile (thin 
lines) of the |A0| distribution as a function of At for the 
small source (a — 0.02i?g, solid) and the large source 
(a = 0.16i?E, dotted). All curves show basically the same 
behaviour: the |A0| values first increase with increasing 
At, then flatten out. This behaviour shows that many or 
most jumps in magnitude and position are dominated by 
one fold-caustic crossing. The time scale is dominated by 
Across- For a slightly larger time interval, the offset does 
not increase significantly any more. Only for much larger 
At's allowing for additional caustic crossings, would an- 
other increase in \A8\ be expected. Offsets of more than 
100 microarcseconds could indeed be reached this way, but 
the characteristic time scale would be depressingly large 
(many decades!). 

The qualitative behaviour of the cases with external 
shear displayed in Figure 1131 is similar to those without 
shear. However, the expected offset values are significantly 
higher here, as already seen in Figs. ITHI and ITT1 again, the 



cases k — 7 = 0.4 and 0.6 appear most promising (sec- 
ond/third row), with median values of between 1 9e and 
2 8e, and 95th-percentiles of 4 8e or higher. These values 
translate into about 7 to 15 microarcseconds (median) and 
> 30 microarcseconds (95th percentile) when applied to 
the quadruple quasar Q2237+0305 (cf. Table [TJ. In fact, 
brightness fluctuations of more than one magnitude have 
been measured in Q2237+0305 on time scales of a few 
months (Wozniak 2000a, b), so these events do occur and 
seem not to be very rare. 

4.3. Discussion and comparison with previous work 

We explored astrometric microlensing effects on short 
timescales (of order months to years) for a set of param- 
eter values k and 7. We find that the relevant time scale 
for measuring relatively large jumps of occasionally many 
tens of microarcseconds can be as short as a few months. 
Such sudden changes of position - produced by caustic 
crossings - are statistically related to fluctuations in the 
apparent brightness of the quasar. Therefore, a good strat- 
egy for detecting this centroid shift would be to measure 
the positions of the most promising lensed quasars very ac- 
curately once, then to monitor the quasars in the optical, 
and - when a photometric microlensing event is detected 
- to perform one or a few more accurate astrometric mea- 
surements. 

We found that the effect is most pronounced for values 
of the surface mass density and shear k — 7 = 0.4 or 
0.6. These parameters happen to be applicable to many 
of the lensed quasar images. The most favorable case is 
the quadruple lens Q2237+0305: because of the closeness 
of the lensing galaxy, the time scale is relatively short (cf. 
also Wozniak 2000a, b) . We also investigated the positional 
shift of the image as a function of source size: whereas a 
quasar with a typical size of 0.16 9e produces median 
offsets of order Be (and 95th-percentiles of about 2 6 E ), 
smaller sources (0.02 6e) reach median values of 2.5 Qe 
and 95th-percentiles larger than 5 0e (where in the case 
of Q2237+0305, 6 E « 7.3 x y/M/M Q hf 5 5 ^arcsec, cf. 
Table HJ. 

Lewis and Ibata (1998) had investigated the astromet- 
ric microlensing effect specifically on Q2237+0305. They 
had found that substantial image shifts of rj 100/iarcsec 
are possible within months. We can confirm this in some 
rare cases. For typical caustic crossings with photomet- 
ric fluctuations of about 0.5 mag, we find values between 
20 and 40 /uarcsec. Salata & Zhdanov (2003) have looked 
into the question of rms fluctuations of the quasar posi- 
tion; however, they only considered large sources (Q.5Re) 
and cases with k + 7 < 0.8. 



4.4. Prospects for observations 

Ground-based differential astrometry in the near in- 
frared is able to achieve measurement uncertainties of 
better than 10 microarcsec, as was reported very re- 
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Fig. 12. The median (thick lines with cross symbols) and 
95th percentile (thin lines with square symbols) of the 
\A8\ distribution as a function of At for |Ato| > 0.2 (left 
panels) and \ Am\ > 0.5 (right panels). The solid and dot- 
ted lines represent the two source sizes a — 0.02Re and 
0.16i?_E, respectively. From top to bottom, the panels cor- 
respond to increasing surface mass density: k = 0.2, 0.4, 
0.6 and 0.8, respectively, with no external shear (7 = 0). 

cently from the Palomar Testbed Interferometer (Lane & 
Muterspaugh 2003). This is only feasible for bright ob- 
jects so far, but it is a very exciting result which opens 
up promising opportunities for the coming years. Another 
instrument promising extremely high astrometric preci- 
sion in the near future (with planet detection as one of its 
main scientific drivers) is the PRIMA instrument at the 
ESO VLTI. It should become efficient in 2004 (Paresce et 
al. 2003; see also Delplancke et al. 2001). The current goal 
at ESO is to achieve 50 /^arcsec accuracy with PRIMA in 
the H and K bands in 2005-2008 and 10 /xarcsec accuracy 
in 2008-2010 (Th. Henning, private communication). 

A number of space-based astrometric projects are 
also underway: The Space Interferometric Mission (SIM) 
is a five year mission scheduled for launch in 2009 
I http://sim.jpl.nasa.gov). The mission's goal is to reach 
an astrometric accuracy of about 1 microarcsecond for a 
predefined grid of objects brighter than 13 mag in the 
visible, which is quasi-inertially tied to a set of distant 
QSOs. For slightly fainter objects (the four images of the 
lensed quasar Q2237+0305 have magnitudes ~ 17), SIM 
is expected to yield 4 microarcsecond absolute positions 
(Unwin et al. 2002). In order to detect center-of-light off- 
sets for multiple quasars, in fact only relative astrome- 
try between the quasar images is required. SIM will also 
be able to measure such positional shifts as a function 
of color, which will give us hints on the physical struc- 



Fig. 13. Same as Fig. El nere for the cases with ex- 
ternal shear 7 — k = 0.2,0.4,0.6,0.8 (from top to bot- 
tom): median (thick lines, crosses) and 95th-percentile 
(thin lines, squares) of |A0| distribution as a function of 
At for \Am\ > 0.2 (left panels) and |Ara| > 0.5 (right 
panels); solid lines for small source (a = 0.02Re), dotted 
lines for large source (0.16-Re), respectively. 

ture of the continuum emission region: presumably the 
cooler/redder part is more extended than the hotter/bluer 
part, which means that we expect larger changes in the 
center-of-light at shorter wavelengths. The GAIA satellite 
is an ESA mission currently scheduled for launch in June 
2010 (see http://sci.esa.int/gaia, Perryman et al. 2001, 
Perryman 2002). With a nominal precision of a few mi- 
croarcseconds for bright objects (about 10 microarcsec- 
onds for 15 th mag objects), it will measure accurate po- 
sitions of 500 000 quasars. So GAIA is expected to be an 
extremely useful instrument for astrometric microlensing 
purposes. It will provide many positional shifts of quasar 
images, along with their lightcurves in many filters. 

5. Summary and Conclusions 

We analyzed the shifts in the center-of-light positions of 
gravitationaly lensed quasar images and the corresponding 
flux variability due to the microlensing effect of stars in 
the lensing galaxy. We found the following results: 

1. The center-of-light of multiply imaged quasars is ex- 
pected to vary on time scales of months to years, due 
to astrometric microlensing effects. 

2. This effect depends on the size of the lensed quasar, 
which we modelled with a Gaussian profile: it is larger 
for smaller sources. This means that measuring cen- 
troid shifts will allow us to constrain the size of the 
quasar. 
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3. We studied eight cases, four without external shear 
(k = 0.2,0.4,0.6,0.8; 7 = 0.0) and four with shear 
( K = 7 = 0.2,0.4,0.6,0.8). The strongest effects are 
seen in the cases k = 7 = 0.4 and 0.6, which in fact are 
very typical parameter values for many lensed quasar 
images. 

4. The effect of center-of-light shift of quasars is not lim- 
ited to multiply imaged quasars: single quasars may 
occasionally be microlensed too (as was suggested by 
Hawkins 1996). The offsets should be relatively small, 
however, due to the fact that the values of K and 7 
will be lower for single quasars than for multiply im- 
aged ones. Our simulation with k = 0.2 comes closest 
to describing these cases with very low convergence. 

5. The centroid shift effect is statistically strongly corre- 
lated with changes in the apparent brightness of the 
quasar. Both occur during the crossing of a (micro- 
caustic: photometric fluctuations > 0.5 mag produce 
offsets > 29e in 50% of the cases, and > 59 e in 5% 
of the cases. These values correspond to about 15 and 
35 microarcseconds, respectively, for the lensed quasar 
Q2237+0305. 

6. The best strategy for measuring centroid shifts is 
therefore to monitor the apparent brightness of multi- 
ple quasars with ground based telescopes and to de- 
termine the relative positions of the various quasar 
images occasionally, ideally before, during and after 
a high magnification event. 

7. In the cases with external shear, the center-of-light 
shift depends strongly on the direction of the shear, 
which can hence be determined in a statistical sense. 

With the next generation of astrometric instruments pro- 
viding an accuracy of order 10 microarcseconds, the astro- 
metric microlensing effect of stars acting on background 
quasars will become detectable. 
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